Time-dependent occupation numbers in reduced-density-matrix functional theory: 
Application to an interacting Landau-Zener model 
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We prove that if the two-body terms in the equation of motion for the one-body reduced density 
matrix are approximated by ground-state functionals, the eigenvalues of the one-body reduced 
density matrix (occupation numbers) remain constant in time. This deficiency is related to the 
inability of such an approximation to account for relative phases in the two-body reduced density 
matrix. We derive an exact differential equation giving the functional dependence of these phases in 
an interacting Landau-Zener model and study their behavior in short- and long-time regimes. The 
phases undergo resonances whenever the occupation numbers approach the boundaries of the interval 
[0, 1] . In the long-time regime, the occupation numbers display correlation- induced oscillations and 
the memory dependence of the functionals assumes a simple form. 
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I. INTRODUCTION 

The effective single-particle Schrodinger equations in 
reduced-density-matrix functional theory [1] (RDMFT) 
differ from the Kohn-Sham (KS) [2] and Hartree-Fock 
(HF) equations in that essentially all orbitals have frac- 
tional (0 < rii < 1) as opposed to integer occupation 
numbers. The freedom to occupy orbitals fractionally is 
helpful in describing strongly correlated systems, where 
strong quantum fluctuations among configurations can 
cause the orbitals to have average occupations that dif- 
fer significantly from or 1. The KS and HF equations 
attempt to reproduce certain observables, the density 
and energy, respectively, with only a single configuration 
(Slater determinant). As a consequence, the orbitals can 
lose any resemblance to the optimal orbitals for describ- 
ing the wave function, the so-called natural orbitals [3]. 
In dynamical problems, time-dependent occupation num- 
bers represent changes in the degree of correlation [4, 5]. 
The linear response of the occupation numbers has been 
shown to be crucial for describing double excitations in 
linear response theory [6, 7]. 

In RDMFT, the wave function is interpreted as a 
functional of the one-body reduced density matrix (one- 
matrix) 7 (1, 1', t) = {^(t)\ft{V)j(l)\*{t)) (1 = r 1 ,a 1 ). 
The equation of motion is (in units \e\ = H = m = c = 1) 

idtl= [\(P~ A) 2 + w,7] + iu, (1) 

where v and A describe time-dependent external electro- 
magnetic fields. In the spatial representation, u is 

(l|fi|l') = j Jd2[vc(l,2)-v c (l',2)} r(12,l'2,t). (2) 

Here, r(12,l'2',t) = |(*(t)|^t(i')^t(2')^( 2 )^(l)|*(t)) 
is the two-body reduced density matrix (two-matrix) and 
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vc is the Coulomb potential. Equation (1) can be closed 
by interpreting the two-matrix as a functional of the one- 
matrix and the initial state. In fact, it follows [8, 9] from 
the Runge- Gross theorem [10], or its extension [11, 12], 
that there exists an exact two- matrix functional r([7], £), 
where we have suppressed the initial state dependence. 
In general F([ 7 ],i) is a memory-dependent functional; 
i.e., it depends on 7 (t') for all t' < t. Instead of propagat- 
ing the one-matrix directly, it may be more convenient 
to propagate its eigenfunctions and eigenvalues, called 
natural orbitals and occupation numbers, respectively, 
according to the equations [4, 8, 13] 

*I0*> = [^(P- Af+v + v ee ] \</> k ), (3) 

nk = (<t>k H <t>k) 

= ilm^TijuVkiij (4) 

ijl 

where the dot represents the time derivative, v ee jk = 
(4>j\v ee \(f)k) = iujk/{n k - rij) (j ^ k), and T ijk i and V ijk i 
are the two-matrix and Coulomb integral expressed in 
the natural orbital basis, respectively. 

In this paper, we report four fundamental results con- 
cerning the time dependence of the occupation numbers: 
(i) a proof that the occupation numbers are time indepen- 
dent when the exact functional r([7],t) is approximated 
by the adiabatic extension of any ground-state (gs) func- 
tional r[ 7 ], (ii) an explicit differential equation for the 
exact memory-dependent functional r([ 7 ],t) in an inter- 
acting generalization of the Landau-Zener (LZ) model, 
(iii) the identification and characterization of correlation- 
induced oscillations in the occupation numbers, and (iv) 
the identification of a universal resonance phenomenon 
responsible for maintaining the Pauli exclusion principle 
in real-time dynamics. Result (i) establishes the need 
for memory-dependent approximations to r([ 7 ],i), while 
(ii-iii) provide exact formulas for memory dependence in 
an important generic case, namely, the crossing of two 
single-particle levels coupled by interactions. 
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II. DEFICIENCY OF THE ADIABATIC 
EXTENSION APPROXIMATION 

If the external potentials v(r, t) and A(r, t) are slowly 
changing functions of time, the wave function remains 
close to the instantaneous ground state. Therefore, it is 
natural to approximate r([7],i) in Eqs. (3) and (4) by 
the gs functional rpy], which we refer to as an adiabatic 
extension approximation (AEA). However, the AEA was 
found to have the deficiency that if a HF-type gs func- 
tional is used, the occupation numbers remain constant 
in time [4, 5, 14]. Since this result relied on the special 
form of HF-type functionals, it left open the possibility 
that time-dependent (td) occupation numbers could be 
obtained with more general gs functionals. Recently, it 
was stated that the occupation numbers are always con- 
stant in the AEA [15], regardless of the gs functional 
that is used. It is important to know whether this state- 
ment is true because it has implications for the design of 
functionals capable of changing the occupation numbers. 
Although the arguments given in Ref. 15 to support the 
statement are incorrect [16], the statement is indeed true. 
We now present a simple proof. 

Consider a system in its ground state at t = to that 
experiences the external driving v(r, t) and A(r, t) for 
t > to. If we exclude external driving that turns on 
discontinuously, then 7(to) = because the system is in 
a stationary state at t = t - This implies nfc(to) = 0. 
Since the AEA is based only on the gs functional T[y], 
at any later time the right-hand side (rhs) of Eq. (4) in 
the AEA is the same as the exact rhs we would have for 
a system just starting in its gs at that time. Thus, it 
vanishes because the gs is a stationary state. We have 
assumed that the j(t) obtained in the AEA remains gs 
ensemble w-representable by a local or nonlocal external 
potential v(ra, rV). If it did not, T[-f] on the rhs of 
Eq. (4) would become ill-defined. 

The natural orbitals are not constant in the AEA since 
they are driven by the external fields. The arguments 
used in the proof do not apply to Eq. (3) because its rhs 
cannot be interpreted as a gs functional. There is a mis- 
match due to the presence of v(r, t) and A(r, t). Namely, 
the gs Y$>) uniquely determined [1] by the instantaneous 
7(t) is generally not the instantaneous gs \ correspond- 
ing to v(r, t) and A(r,t). 

In contrast, the occupation numbers are not driven 
directly by the external fields because v(r,t) and A(r, t) 
do not appear in Eq. (4) . Instead, they are driven purely 
by the internal correlation of the system, as only the 
correlation part of the two-matrix gives a nonvanishing 
contribution on the rhs of Eq. (4) [4, 5]. The correlation 
part of the two- matrix is defined as T c = T — Thf, where 
T HF (11', 22') = 7 (1, l') 7 (2, 2')-7(l, 2')7(2, 1') is the HF 
two-matrix. On the basis of the above proof, we can make 
the stronger statement that the are driven purely by 
nonadiabatic correlation, i.e., the difference between T c 
and its instantaneous gs value. 

What is T[7] missing that makes it incapable of gener- 



ating hk 7^ in Eq. (4)? For two-electron singlet states, 
the only difference between the exact functional r([7],t) 
and r[7] are relative phases that correspond to the rel- 
ative phases between the configurations that compose 
the wave function [13]. These relative two- matrix phases 
must differ from their gs values to yield hk ^ in Eq. (4) . 
The AEA fails because, being based solely on the gs func- 
tional r[7], it cannot change the two- matrix phases away 
from their gs values. In the general TV-electron case, the 
functional r([7],t) differs from r[7] in more degrees of 
freedom than just these relative phases. An alternative 
functional theory that might allow one to take into ac- 
count the relative two-matrix phases in an effective way 
has been introduced [15]. The phases of the natural or- 
bitals, which in RDMFT are undefined, are incorporated 
into the basic independent variable of the functional. 

The relative phases can be seen explicitly by consider- 
ing the Lowdin-Shull [17] wave function for two-electron 
singlet states, which can be written as 

l* LS ) = ^e-'^e-^V^a^l), (5) 

where a\ a are the natural orbital creation operators and 
Cfc are the relative phases between the configurations. 
The corresponding expression for the two-matrix in the 
basis of natural orbitals, T^i — \e l2 ^ k ~^ \fh~jXik~5ij5ki, 
shows that varying the (k changes the phases of the el- 
ements Tijki, i.e., it changes relative phases in the two- 
matrix. There is a one-to-one correspondence between 
the relative two- matrix phases and the Cfe- The (k are 
not invariant under redefinitions of the natural orbital 
phases, so they can only be uniquely defined with respect 
to a given choice of td natural orbital phases. However, 
for any choice of td natural orbital phases, the Cfc have 
unique gs values. To generate td occupation numbers in 
Eq. (4), the (fc in the functional r([7],t) must differ from 
their gs values. Memory dependence enters r([7],t) ex- 
clusively through the £fc, which are functionals of "f(t). 
The relative phases can be given a geometric interpreta- 
tion, which will be discussed in Sec. III. 

Despite the failure of the AEA in Eq. (4) , it is possible 
to obtain td occupation numbers in an approach based 
only on the gs functional T[y]. In the adiabatic regime, td 
occupation numbers can be obtained on the fly during the 
propagation of Eq. (3) from a condition of instantaneous 
occupation number relaxation (IONR) [13]. The IONR 
approximation reproduces an adiabatic approximation in 
linear response theory [14] , but it is more general because 
it applies to fully nonlinear real-time dynamics. It cap- 
tures the lowest-order nonadiabatic effects even though 
it lacks memory dependence. However, if the occupa- 
tion numbers deviate greatly from their instantaneous gs 
values, the IONR approximation breaks down. In such 
cases, it is necessary to propagate Eq. (4). Recently, an 
approach was proposed in which an approximation for T c 
is obtained from semiclassical propagation of the iV-body 
density matrix. This T c is then used in the propagation 
of Eq. 1, yielding td occupation numbers. 
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III. 



INTERACTING LANDAU-ZENER MODEL 



Consider two generic single-particle states whose bare 
energies cross linearly in time. If we occupy the system 
with two electrons and allow interactions between them, 
we obtain an interacting generalization of the Landau- 
Zcner (LZ) model. We consider the dynamics in the 
sector of spin-singlet states. The extension to the full 
Hilbert space is an interesting problem for future study. 
The spin-singlet sector has three states. Having at least 
three states is essential for representing genuine inter- 
actions. The three levels undergo a correlated avoided 
crossing (see Fig. 4 in Ref. 13), which is one of a hierar- 
chy of multiplet-like avoided crossings that one encoun- 
ters in many-body systems. Occupation numbers vary 
most rapidly near such avoided crossings. The Hamilto- 
nian with the most general one-body and two-body terms 
in the spin-singlet sector is 



H = ^V -ff+U + W, 




(6) 



FIG. 1: Time dependence of A and C for Ui = Ui = 3/2, 
\W\ = 0, and driving potential V = (-2,0, At). 



where a = £ CT (C c\ a )a 
potential, and 



C2a 



V acts as an external 



W ={W X - iW 2 )c\ t c\ i C2iC2 t + {Wi + iW 2 )cl t cl l c n c lf . 

The spin-summed one-matrix in this model is simply 
a Hermitian 2x2 matrix 7 = J + 7 • a, where 7 = 
A(sin 9 cos ip, sin# sin 95, cos 9) is similar to a pseudospin 
vector with Bloch sphere angles 9 and tp, but its modulus 
A is less than 1 due to correlations. The natural orbitals 
are chosen to be (f> a = (cos(<9/2)e-^/ 2 , sin(6>/2)e^/ 2 ) T 
and <j) b = (-sin(6>/2)e-^/ 2 ,cos(6>/2)e^/ 2 ) T . The cor- 
responding occupation numbers are n a = 1 + A and 
rib = 1 — A. Equations (3) and (4) become 



i<t>k = {^V ■ a + v ee )(f> k 



A = -(u ao 



Ubb)- 



(7) 

(8) 



In Eq. (7), v ee is the contribution of U and W to the 
effective single-particle Hamiltonian. In the natural or- 
bital basis, its off-diagonal elements are v ee ,ab = vl e ha — 
—iu a b/2A, while its diagonal elements are indeterminate 
because the phases of the natural orbitals are undefined. 
Direct calculation with |5 ,LS ) in Eq. (5) gives 

Uab = -iU{\ + Ber ac - ) sin 9 cos 9- iAUA sin 9 
+ i\W\(l + Be- t2C ) sin 9 cos 9 cos(2<^ - to) 
-\W\{1- Be- i2C )sm9sin(2ip-uj), (9) 



where B = VI - A 2 , \W\e~ 1 " = W x - 
U 2 )/2, AU = (C/i - U 2 )/2 and ( = ( a 



iW 2 , U =(U 1 + 
- (b is the relative 



two-matrix phase. The diagonal elements of u enter in 
Eq. (8). We have u aa = —Ubb and 

u aa = -UB sin 2 9 sin 2( - 2\W\B cos 6>sin(2^-u;) 
x cos 2( - I W\B(1 + cos 2 9) cos(2ip - u) sin 2C,. 

The dynamical equation for £ can be derived from the 
stationarity of the action S[&] = dt(^\H - id t \^} 
with respect to A. We find 

C = \v -7 + ^4 sin26,cos2 C + AUcos9 + 

A 2 B 2 B 

x (l + cos 2 6>)cos(2^-w)cos2C- |V^|— cos 6» 

B 

x sin(2<^ — to) sin 2( — ip cos 9. (10) 
Here, we have used = V ■ 7 + + W[*] with 

U[V} = ^(1 + cos 2 9) - jBsin 2 9cos2( + AUAcos9, 
W[^\ = ^ sin 2 9 cos(2<^ - u) + \ W\B cos 9 sm(2p - u) 
x sin 2C - ^y^ 1 + cos2 e ) cos(2p - to) cos 2( 

and i dt(^\^) = dt(fi + Aip cos 6 + AC,). 

As Eqs. (7), (8) and (10) describe all of the degrees of 
freedom of \^> LS ) in Eq. (5), they are equivalent to the 
Schrodinger equation. In RDMFT, Q is interpreted as 
a functional of 7(i). If the mapping V(i) — > j(t) is in- 
vertible for a given initial state, the memory-dependent 
functional C([^f],t) is given uniquely by the solution of 
Eq. (10). In Fig. 1, A and ( are plotted for a linear- 
time driving potential V = (— 2,0, At). The initial state 
is the instantaneous ground state at to = — 4. Although 
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in numerical simulations it is not possible to specify the 
initial condition at t = — oo as in the LZ model, we 
chose an initial time early enough so that the td A and 
C are numerically converged with respect to the limit 
to — » — oo. In the limit t — > oo, A displays persistent 
oscillations with frequency 2U. This is surprising be- 
cause in the long-time limit, where the system is a non- 
stationary state, one expects the oscillations to have the 
frequencies Qij = Ei — Ej, where Ei are the adiabatic 
energy levels. The components of 7 do indeed oscillate 
with the frequencies fiy. Due to the divergence of the 
linear-time driving potential V, the f2,j diverge in the 
limit t -> 00 (E-l ~ -V 3 + U, E 2 ~ and E 3 ~ V 3 + U). 
In contrast, A, which describes the occupation number 
degrees of freedom, oscillates with the constant frequency 
^32 — ^21 = 2f7. This expression demonstrates that 
all three states are participating in the oscillations; they 
have no analog in two-state systems. The emergence of 
the Hubbard energy U is a consequence of the fact that 
the occupation numbers arc driven exclusively by the in- 
ternal correlation of the system. One can expect to see 
correlation-induced oscillations in n k whenever a many- 
body system traverses a correlated avoided crossing. 

We also observe that ( decreases linearly in the limit 
t — > 00. The asymptotic slope, —U, is exactly half the 
frequency of the oscillations in A. The relationship be- 
tween A and £ in the long-time regime will be discussed 
further in Sec. Ill A. Superimposed on the linear depen- 
dence are nonlinear oscillations with a "sawtooth" pat- 
tern. The rapid jumps in £ come from the term tpcosO 
in Eq. (10). To understand this, consider the dynam- 
ics of 7 in three-dimensional space. The vector 7 be- 
gins at the north pole (0 = 0) of the Bloch sphere at 
t = —00. Then, it follows adiabatically the driving vec- 
tor — V as it rotates toward the south pole. Although 
—V/\V\ approaches the south pole in the limit t — > 00, 7 
does not. Instead, due to nonadiabatic transitions, it pre- 
cesses continuously around the south pole with constant 
73 iOC = 73(00). However, before it reaches this asymp- 
totic behavior, there is an interval of time during which, 
periodically, 7 passes close to the south pole (see Fig. 1 
in Ref. 13). For each such time, ip is strongly peaked, 
thus inducing jumps in (. 

The phases can be given a geometric interpretation. 
Consider a cyclic evolution on the time interval [0,T]. 
A cyclic evolution is one for which \^(T)) differs from 
|*(0)) by only an overall phase, i.e. |*(T)> = e _i,/ |*(0)). 
The geometric phase [18] for the general two-electron 
spin-singlet state in Eq. 5 is 

ij dt{£\i) = ^<j>{in k (<p k \d<pk) + n k d( k ), (11) 

where = \^ (t)) and fi(t) is any real-valued 

function for which /i(T) — /z(0) = v. The first term in 
Eq. 11 is the geometric phase associated with the nat- 
ural orbitals. The second term shows that the time- 
evolving Cfe give a contribution to the total geometric 
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FIG. 2: Example of a resonance in £ when nt approaches 1. 
Parameters are Ui = U 2 = 1, \W\ = 0, and V = (-2, 0, 2t). 



phase above and beyond the natural orbital contribution. 
In the present model, the geometric phase simplifies to 
§(Acos6dtp + AdQ. In the noninteracting case, this be- 
comes § cos 6dtp, which is just the familiar result for the 
pseudospin of a two-level system. Interactions modify 
the noninteracting geometric phase in two ways. First, 
the factor A, which is between and 1, reduces the nat- 
ural orbital contribution. From its definition, we have 
A = (n a — n b )/2, where n k = n k ^ + n ki and n a > n b . 
A = 1 in the noninteracting case, and it decreases as 
the correlation of the state increases. Second, interac- 
tions introduce the additional term § n k d(, k = § AdQ. 
This term vanishes in the noninteracting case because for 
n k = const, § n k dQ k = n k § d( k = 0. Interactions have a 
similar effect on the gs Berry phase. 

It is important to remember that Q is a relative phase 
that depends on our choice of td natural orbital phases. 
Redefining the natural orbital phases so as to satisfy the 
parallel transport condition Im((f) k \<f) k ) = redefines ( 
according to £ — > £ — ip cos 9. Then, the geometric phase 
in Eq. (11) becomes i dt(£\£) =Y^k§ n kd( k = § AdC,. 

There is a remarkable aspect of the ( k worth discussing. 
When one of the n k approaches one of the boundaries of 
the interval [0,1], the phase (,k undergoes a resonance, 
jumping by 7r/2. The effect of this resonance is to change 
the sign of h k , thereby keeping the occupation numbers 
in the allowed interval [0, 1]. The upper bound is a con- 
sequence of the Pauli exclusion principle. 

The width of the resonance depends on how close n k 
comes to the boundary. Figure 2 shows A and £ in an 
interval where A approaches close to 1, i.e. n a ^ = n a ^ — > 
1 and nff = nbi — > 0. We can obtain a universal equation 
for £ near the resonance by keeping only the leading terms 
in the limit A — > 1, which is equivalent to B — > 0. From 
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Eqs. (8) and (10), we find 



o. 



(12) 



Integration gives £ = const /B 2 , and further integration 
gives the memory-dependent functional C([7],i). Choos- 
ing the case \W\ = for simplicity, the explicit solution 
of Eq. (12) is 



C= ^tan-^t-tc)]. 



(13) 



This is the formula for a resonance of width or 1 centered 
at t = t c . We expect a similar resonance phenomenon 
will occur in TV-electron systems at the boundary rik = 1. 
The boundary n^ — might be more subtle, as it is an 
accumulation point of the spectrum of 7. 



A. Long-time regime 

We consider the long-time limit of the model when it 
starts in the ground state at t = —00 and experiences 
the external driving V = (Vi,0,i/r). For simplicity, we 
set U = U\ = U2 and \W\ — 0. The exact asymptotic 
behavior of A(t) in the limit t — >■ 00 is 



A(t) = \jA 2 + A 2 cos [2Ut - (632 - e 2 i ) 



(14) 



where Vee.oo = lim t ^oo V^^] acts as an effective Hamil- 
tonian and A and £ appear as canonically conjugate 
action-angle variables. Equations (17) and (18) are inte- 
grable and V ee .oo is a constant of the motion. Since V ee>oc 
can be expressed as a function of A, £ and 73,00, there is 
the following instantaneous relationship between A and 
( in the long-time limit: 



cos2£ 



2A 2 (l - 



04 2 -7f,oo) 



y/l - A 2 (A 2 



(19) 



Therefore, the 7 dependence of C([j},t) separates into 
two distinct types: (i) an ultra-local (instantaneous) de- 
pendence on A(t) and (ii) an ultra-nonlocal dependence 
on 7(f) near t = that enters only through the con- 
stant 73,00- The dependence on 73,00 contains informa- 
tion about the nonadiabatic transitions that occurred 
near t = 0. In the long-time regime, it can be viewed as 
initial-state dependence specified after the avoided cross- 
ings are complete. Memory dependence becomes simple 
because nonadiabatic transitions are no longer occurring 
in the long-time regime. Equations (15-18) suggest that 
memory dependence among certain degrees of freedom 
will be weaker when their mutual dynamics are nearly 
integrable and nearly decoupled from other degrees of 
freedom. 



B. Short-time regime 



where A 2 



73,oo 



+2|c 2 | 2 (| Cl | 2 + | C 3| 2 ),A 2 ^4| Cl 



M 2 M 



@ij = Arg(ci/cj) and Cj are the coefficients in the expan- 
sion I*) = J2Li QCxpH / dt(Ei - i<$i|at$i»]|*i) over 
the adiabatic eigenstates Since the pseudomagnetic 
field V diverges as t — > 00, the pseudospin-like vector 7 
precesses more and more rapidly around the south pole 
of the Bloch sphere. As t — > 00, the azimuthal angle 
ip grows quadratically in time, i.e., ip ~ t 2 /2r. Hence, 
any terms in the dynamical equations containing peri- 
odic functions of ip quickly average to zero. Dropping 
such terms in Eqs. (8) and (10), we find 



A = -[LB sin 2 (9 sin 2C 

c 



l + Bcos2C 2 

U ; COS - 



UA 
~2B 



(15) 

sin 2 9 cos 2(. (16) 



Although 9 appears in these equations, it can be elim- 
inated in favor of A by means of the asymptotic rela- 
tionship 73,00 = Acos9. Therefore, asymptotically the 
equations for A and £ decouple from those for the or- 
bital variables 9 and <p. Equations (15) and (16) can be 
expressed in the form of Hamilton's canonical equations, 



dV e( 



C = 



,00 



dA 



(17) 
(18) 



We now study £ shortly after an external driving po- 
tential is turned on for a system in the ground state at 
finite t = to. The driving potential is taken to have the 
form V = (Vi, V2, Vs(t)), where V\ and V2 are constant 
and Vz(t) is an arbitrary continuous function for which 
Vz{to) 7^ 0. This is analogous to a td local external po- 
tential in a continuous system. The stationary conditions 
imply 7 (to) = and C(*o) = 0. The lowest nonvanish- 
ing time derivative at the initial time is (pita) = V-j^o)- 
This induces nonvanishing third time derivatives for 9, 
(p, and £. The lowest nonvanishing time derivative of A 
is of fourth order. The relative phase and the occupa- 
tion numbers change even more slowly if |W| = 0. In 
this case, the lowest nonvanishing time derivative of A 
is the fifth time derivative, and the changes proceed via 
V 3 -> (p -> -> C (4) A^ ] . This means that the exter- 
nal driving induces first a current, the current leads to 
changes in the density, the changes in the density cause 
changes in (, which, finally, induce changes in the occu- 
pation numbers. 



IV. CONCLUSIONS 

The occupation numbers are important degrees of free- 
dom that provide information about the correlation of 
a many-body state. Their dynamics are determined by 
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the equation of motion for the one-body reduced den- 
sity matrix. This equation of motion contains the two- 
body reduced density matrix T(t), and we have proved 
that approximating T(t) by the adiabatic extension of 
the ground-state functional Fby] does not generate time- 
dependent occupation numbers. For two-electron sys- 
tems, this deficiency can be explicitly traced to the 
fact that the ground-state T[j] lacks dynamical relative 
phases Qk present in the exact T(t). Physically, the 
describe nonadiabatic interaction effects. The variables 
rik and (k are canonically conjugate variables that give 
a contribution to the geometric phase above and beyond 
the contribution from the parallel transport of the natu- 
ral orbitals. Additionally, they have the important func- 
tion of maintaining compliance with the Pauli exclusion 
principle. 

We derived an exact differential equation for the £fc 
in a generalization of the Landau-Zener model that in- 
cludes interactions. Introducing two-body interactions 
splits the single Landau-Zener avoided crossing of two 
levels into a set of three pairwise avoided crossings among 
three levels. After the system traverses these correlated 
avoided crossings, the occupation numbers display oscil- 
lations whose frequency is determined by the Hubbard 
interaction U . These correlation-induced oscillations, 
which depend on the participation of all three states in 
an essential way, are an observable effect of interactions 
in quantum many-body dynamics. 

How to account for memory dependence in time- 
dependent density functional theories, especially in real 
time, is an important unsolved problem. Equation (10) 
gives an explicit example of a differential equation whose 
solution determines the memory-dependent functional 
C([l],t) and thereby the complete functional r([7],£). 
We have found that the 7 dependence assumes a sim- 
ple form in the long-time regime, comprising both an 
instantaneous dependence and an ultra-nonlocal depen- 
dence representing past nonadiabatic transitions. Mem- 
ory dependence also simplifies within the so-called in- 
dependent crossing approximation, where each avoided 
crossing is treated as independent of the others. At this 
level of approximation, memory dependence enters via 
dynamic and geometric phases and the amplitudes a of 



the adiabatic eigenstates [13]; the which are assumed 
to be constant for all states except those participating in 
a given avoided crossing, retain the memory of past nona- 
diabatic transitions. Conversely, the exact memory de- 
pendence is nontrivial while nonadiabatic transitions are 
occurring, e.g., while the system is traversing an avoided 
crossing. 

Two adiabatic energy levels that undergo an avoided 
crossing in real time generally intersect at certain points 
in the complex plane of time. The functional dependence 
of r([7],i) might simplify at such times, and indeed, the 
contour integral in the Dykhne formula suggests that in 
the adiabatic regime memory dependence might be easier 
to handle in complex time. We also note that Eq. (10), 
viewed as a differential equation in the complex plane 
of time, has singularities when n a — njj. In the inter- 
acting Landau-Zener model, we have found that the in- 
stantaneous ground-state n a and n?, have intersections in 
the neighborhood of the intersections of E\ and Ei [19]. 
Hence, the occupation numbers contain useful informa- 
tion about the avoided crossings of the adiabatic energy 
levels. 

The interacting Landau-Zener model studied here 
might provide insight into the nonadiabatic dynamics of 
many-body systems. In situations where a pair of nat- 
ural orbitals are strongly coupled and the others can be 
treated as an approximately inert background, it is pos- 
sible to map the subspace dynamics onto the interacting 
Landau-Zener model. The model parameters V, U and 
W, representing the effective "screening" of the other 
degrees of freedom, will be time dependent. It might be 
possible to devise a memory-dependent functional based 
on Eq. (10) and the assumption that such pairwise in- 
teractions of the natural orbitals can be treated one at a 
time. 
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